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;_( ■ Abstract 

■^ . We present analytic representations of solutions to the 2-photon and two-mode 

C^ . quantum Rabi models in Bargmann spaces of entire functions. For each case, we 

obtain a transcendental function whose zeros give the energy spectrum of the model. 

The zeros can be found numerically by standard root-search techniques. 
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1 Introduction 

Spin-boson systems describe interactions between spins and harmonic oscillators (boson 
^ ! modes) and have played a prominent role in modeling the ubiquitous matter-light inter- 

t^^ . actions in modern physics. One of the well-known spin-boson systems is the quantum 

"^ ■ Rabi model, introduced by Rabi over 70 years ago. This model describes the interaction 

cn ' of a two-level atom with a harmonic mode of quantized electromagnetic field. Due to the 

simplicity of its Hamiltonian, the Rabi model has served as the basis for understanding 
matter-light interactions, and has a variety of applications ranging from quantum optics 
[1] to solid state semiconductor systems [2] and molecular physics [3]. It has also played a 
significant role in the novel research field of cavity and circuit quantum electrodynamics 
[4, 5]. 

The Hamiltonian of the Rabi model is given by 

HR = ub^b + Aa, + ga^ib^ + b), (1) 

where g is the spin-boson interaction strength, cr^, (Tx are the Pauli matrices describing the 
two atomic levels separated by energy difference 2A, and 6^ (6) are creation (annihilation) 
operators of a boson mode with frequency u. In the Bargmann representation of the boson 
creation and annihilation operators, in terms of the two-component wavefunction ip{z) = 
{'ip+{z),ilj-{z)) , the time-independent Schrodinger equation of the model gives a system 
of two coupled linear differential equations with rational coefficients. The system depends 



on a spectral parameter E. The wavefunction components are elements of Bargmann- 
Hilbert space of entire functions of one variable z G C [6]. The scalar product of any two 
elements f{z),g{z) in the Bargmann-Hilbert space is given by 



{f,9)=J f{z)g{z)dfi{z), (2) 

where djj,{z) = -e"'^' dxdy. The energy E belongs to the spectrum of the problem if 
and only if for this value of E the system of coupled differential equations (given by the 
time-independent Schrodinger equation) has entire solution [6]. 

By means of the Bargmann-Hilbert space approach, Braak [7] recently presented a pair 
of transcendental functions defined as infinite power series expansions with coefficients 
satisfying three-term recurrence relations, and argued that the spectrum of the Rabi 
model is given by the zeros of the transcendental functions. This theoretical progress has 
renewed the interest in the Rabi and related models [8, 9, 10, 11, 12, 13]. 

In this work we analyze the 2-photon and two-mode generalizations of the quantum 
Rabi model within the framework of Bargmann-Hilbert spaces of entire functions. The 
2-photon and two-mode Rabi models are phenomenological models describing a two-level 
atom interacting with 2 photons and 2 harmonic modes, respectively. They can be exper- 
imentally realized in circuit quantum electromagnetic systems [5], and have established 
applications in many research fields, including Rubidium atoms [14] and quantum dots 
[15, 16]. We present analytic representations of solutions to the two models. For each 
case, we analytically determine a transcendental function whose zeros give the energy 
spectrum of the model. In the Appendix, we revisit the driven Rabi model with broken 
Z2 symmetry and derive a continued fraction equation for its energy spectrum. 

2 2-photon quantum Rabi model 

The Hamiltonian of the 2-photon Rabi model reads 

H2p = cob^b + Aa^ + g a^ [{b^f + 6^] . (3) 

Let us make a canonical Bogoliubov transformation from 6, 6^ to the squeezed bosons 

a, a\ 

a-\- T a^ ^^ T a + a^ 

b = , b^ = , (4) 

where |r| < 1 is a real parameter. In terms of the squeezed bosons, the Hamiltonian (3) 
takes the form 



H2p = Aaz + ^ 

1 — r^ 



cur + ga^{l + T^)) ({a^f + a^ 



1 

+ {uj{l + T^) + Agra^ ) a} a + tor^ + 2gTa^ . (5) 



Introduce the operators K±^ Kq 

1 



K, 



-[a 
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t^2 



2 ' 



K. 



— a'a H — 
2 V 2 



(6) 



Then (5) becomes 

H2p 



Act, + 



l-r2 



2(cuT + ga,{l + T^)){K+ + K. 



2 Uj(l + T^ 



Agra^ Kq 



-r- 



(7) 



The operators K±,Ko form the su{l,l) Lie algebra. Its quadratic Casimir, C = 
K^K_ — Kq{Kq — 1), takes the particular values C = ^ in the representation (6). This is 
the well-known infinite-dimensional unitary irreducible representation T>^{q) of s-u(l, 1) 
with q = |, |. Thus the Fock-Hilbert space decomposes into the direct sum of two 
subspaces "H^ labeled by g = 1/4, 3/4. 

In the Bargmann-Hilbert space of analytic functions labeled by Bargmann index q, 
the operators K±, Kq (6) have the single-variable 2nd differential realization [17] 



d ^^ z ^^ ^ dp , d ,^, 

A-. = .- + ,. A-, = -. A-_ = 2.-^ + 4,-. (8) 

Using this differential realization, working in a representation defined by a^ diagonal and 
choosing r such that ujt + g{l + r^) = 0, i.e. 






i-Q), n = Ji 






(9) 



where 



2s 



< 1, then the transformed 2-photon Rabi Hamiltonian becomes the matrix 



differential operator 



( 2c^^ H + ^) - 1^ 



H. 



2p 



A 



A \ 

'^^ +l[oji2-n')z-8gq]i 



V 



n^ dz2 ^ n 

_22. 

n 



dz 






(10) 
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In terms of two-component wavefuntion '^/'(-z) = (■?/'+ (z), ■?/'_ (2;)) , the time-independent 
Schrodinger equation, H2pip{z) = EiIj{z), yields a system of coupled differential equations. 



2ujn iz— + q] CO- E 

\ dz 2 



i)+ + Av-. 



0, 



'11^ 



^gzj^ + (-2a;(2 - Q^)z + IQgq) ^ 

+2gz - 2a;(2 - n'^)q + {-u + E)VL 



^_ - VlA ^_, 



0. 



(12) 



This is a system of differential equations of Fuchsian type. Solutions to these equations 
must be analytic in the whole complex plane if E belongs to the spectrum of H2p- So we 
are seeking solutions of the form 

oo oo 

i,^{z) = y: nm z\ i,_{z) = y: k{e) z\ (13) 

n=0 ra=0 

which converge in the entire complex plane, i.e. solutions which are entire. 
Substituting (13) into (11), we obtain 

^^ " E + \uj-{2n + 2q)ujn^'' ' ^^^^ 

So /C^ is not analytic in E but has simple poles at 

E = — uj + {2n + 2q)ujVt, n = 0, 1,---. (15) 

The energies (15) appear for special values of model parameters [18, 19] and correspond 
to the exceptional solutions of the 2-photon Rabi model. If (15) is satisfied, the infinite 
series expansions (13) truncate and reduce to polynomials in z but only if the system 
parameters satisfy certain constraints [19]. Majority part of the spectrum of the 2-photon 
Rabi model is regular for which (15) is not satisfied. The regular spectrum of the model 
is given by the zeros of the transcendental function F{E) obtained below. Thus similar to 
the Rabi case, the spectrum of the 2-photon Rabi model consists of two parts, the regular 
and the exceptional spectrum. 

From (12), we obtain the 3-step recurrence relation for /C~, 

/C-+1 + An /C- + Bn IC-_, = 0, n > 1, (16) 



where 



A-n 



ig{n+ l){n + 2q) 



i2n + 2q)u{2 - n 



2^ 



1 A2 \ 

\ 2 E + \uj - {2n + 2q)uQ J 

^" " 4(n + l)(n + 2g)- ^^^^ 

The coefficients An, Bn have the behavior as n — )■ oo 

/l„~an", 5„~6n^ (18) 

with 

a = -^(2-fi'), a = -l, 6 = i, /3 = -2. (19) 



Thus the asymptotic structure of solutions to the 2nd equation of (16) depends on the 
Newton-Puiseux diagram formed with the points Po{0, 0), Pi(l, —1), -P2(2, —2) [20]. Let 7 
be the slope of PqPi and 6 the slope of P1P2 so that 7 = a and 6 = f3 — a. Then we have 
'-f = 6 = a. The characteristic equation of the n > 1 part of (16), t^ = 0, has two equal 
solutions ti = t2 = 0. Applying the Perron- Kreuser theorem (i.e. Theorem 2.3 of [20]), 
we have 

Jim sup(|/C-|n!)"=0 (20) 

for all non-trivial solutions of (16). The zero limit means that the n > 1 part (i.e. the 
truly 3-term part) of (16) possesses a minimal solution /C™" and the corresponding infinite 
power series expansions ■?/'™"(2;), obtained by substituting /C™" for the /C~'s in (14) and 
(13), converge in the whole complex plane, i.e. they are entire. 

By the Pincherle theorem (i.e. Theorem 1.1 of [20]), the ratio of successive elements 
of the minimal solution sequence /C™"^ is expressible as continued fractions, 

V"min D D D 

p _ ^n+1 _ -°"+l -°"+2 -P«+3 /^-.x 

l^n ^n+1— ^n+2— ^n+3^ 



(22) 



which for n = gives 






R — ^ — 

irmin 


A,- 


B2 B-j, 

A2- A3- 


Note that the ratio Rq = ^^ involves 


irmin 


, although 



although the above continued fraction 
expression is obtained from the 2nd equation of (16), i.e the recurrence (16) for n > 
1. However, for single-ended sequences such as those appearing in the infinite series 

trrnin 

expansions (13), the ratio Rq = t^htt of the first two terms of a minimal solution is 
unambiguously fixed by the first equation of the recurrence (16), namely. 



Ro — —Aq 



I6gq 



2gw(2 -V?)-\E^-uj . ^ 



(23) 



In general, the Rq computed from the continued fraction (22) can not be the same as 
that from (23) for arbitrary values of recurrence coefficients An and -B„. As a result, 
general solutions to the recurrence (16) are dominant and are usually generated by simple 
forward recursion from a given value of ICq . Physical meaningful solutions are those that 
are entire in the Bargmann-Hilbert space [6]. They can be obtained if E can be adjusted 
so that equations (22) and (23) are both satisfied. Then the resulting solution sequence 
/C^(-E) will be purely minimal and the power series expansion (13) will converge in the 
whole complex plane. 

Therefore, if we define the transcendental function F{E) = Rq + Aq with Rq given 
by the continued fraction in (22), then the zeros of F{E) correspond to the points in 
the parameter space where the condition (23) is satisfied. In other words, F{E) = 
is the eigenvalue equation of the 2-photon Rabi model, which may be solved for E by 



standard nonlinear root-search techniques (see e.g. [21, 22] and references therein). Only 
for the denumerable infinite values of E which are the roots of F{E) = 0, do we get entire 
solutions of the differential equations (11) and (12). 

3 Two-mode quantum Rabi model 

The Hamiltonian of the two-mode Rabi model is 

H2m = co{blh + hlh) + Aa, + g(j^{h{h\ + 6162), (24) 

where we assume that the boson modes are degenerate with the same frequency u. In- 
troduce the two-mode Bogoliubov transformation, 

_ai + aa\ ,t 0-02 + 0.1 , a2 + (ja\ 4. aai + a\ , . 

^^"Tr^F' ^^"7r^^' ^^"Tr^F' ^^"Tr^^' ^ ^ 

Here |o"| < 1 is a real parameter and 01,02,01,02 are squeezed bosons satisfying the 
canonical commutation relations [oj,oJ] = 1, [oj,Oj] = [oj,o^] = [oJ,o^] = 0, i,j = 1,2. 
In terms of the 2-mode squeezed bosons, the Hamiltonian (24) has the form 



H. 



2m 



2ua + ga^il + c^) ) {alal + O1O2) 






+ [uj{l + a^) + 2ga a^j (o{oi + 0^02) +2uja^ + 2ga cr^;] + Aa^. (26) 
Introduce the operators K±, Kq 



K^ = o|o2, K_ = O1O2, Kq = -(ojoi + O2O2 + 1). (27) 



Then (26) becomes 



Him = Aa 



2cua + ga^{l + a^)) {K+ + K^ 



1-^2 
+2 (w(l + (T^) + 2gaa^) Kq] - u. (2^ 



The operators K±,Ko form the s-u(l, 1) Lie algebra. Its quadratic Casimir, C = 
K^K_ — Ko{Ko — 1), takes the particular values C = /t(l — n) in the representation (27), 
where k = 1/2, 1, 3/2, ■ ■ ■. This is the well-known infinite-dimensional unitary irreducible 
representation of su{l, 1) known as the positive discrete series V~^{k,). Thus the Fock- 
Hilbert space decomposes into the direct sum of infinite subspaces 'H'^ labeled hj k, = 
1/2,1, 3/2, ■■•. 

In the Bargmann-Hilbert space of entire functions labeled by Bargmann index k, the 
operators K±, Kq (27) have single-variable differential realization [17] 

Ko = z— + K, K+ = z, K_ = z—^ + 2k—, (29) 
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where k = 1/2, 1, 3/2, ■ ■ ■. Using this differential reahzation, working in a representation 
defined by a^ diagonal and choosing a such that 2uj(7 + g{l + o"^) = 0, i.e. 



oo 



a 



1-A), A 



'1- 



u' 



(30) 



where 



operator 



< 1, then the transformed Hamiltonian (28) becomes a matrix differential 



H' 



2m 



I 2ujK(zJ-^+ K] -u 



\ 



A 
-\-& + -A [^(2 - A^)- - 2^-] i 



(31) 



/ 



In terms of two-component wavefuntion (f){z) = {(f)^{z),(f)^{z)) , the time-independent 
Schrodinger equation, H2m(p{z) = E(f){z), yields a system of coupled differential equations, 



2a;A i z—- + k 

\ dZ y 

rf2 



uj-E 



+ + A0_ = 0, 
d 



2c/^-^ + (-2w(2 - A2)z + AgK^ 

+2gz - 2a;(2 - A^)^ + {E + u)A] 0_ - AA 



0. 



(32) 



(33) 



This is a system of differential equations of Fuchsian type. Solutions to these equations 
must be analytic in the whole complex plane if E belongs to the spectrum of H2m- Similar 
to the 2-photon Rabi case, we seek solutions of the form 



n=0 

which converge in the entire complex plane. 
Substituting (34) into (32), we obtain 



E QniE)^" 
n=0 



(34) 






A 



■Q". 



E + UJ - (2n + 2/t)a;A' 
So Q^ is not analytic in E but has simple poles at 

E = -uj + {2n + 2K)coA, n = 0, 1, 



(35) 



(36) 



The energies (36) appear for special values of model parameters [19] and correspond to 
the exceptional solutions of the two-mode Rabi model. If (36) is satisfied, the infinite 
series expansions (34) truncate and reduce to polynomials in z but only if the model 
parameters obey certain constraints [17]. Majority part of the spectrum of the two-mode 
Rabi model is regular spectrum which does not have the form (36). The regular spectrum 
of the model is given by the zeros of the transcendental function G{E) obtained below. 



Thus again, the spectrum of the two-mode Rabi model consists of two parts, the regular 
and the exceptional spectrum. 

From (33), we obtain the 3-step recurrence relation for Q~, 

Qi + Co Qo" = 0, 

Q-+i + CnQ- + DnQ-_, = Q, n>l, (37) 

where 

Cn = —, ^7 [-(2n + 2/s:)w(2-A2) 

2^(n+l)(n + 2/t) L ^ ' ^ ' 

I A2 \ ^ 

^ V ^^~ E + oo- {2n + 2K)ujk) 

Dn = -7 -. (38) 

(n + l)(n + 2/t) ^ ' 

The coefficients C„, D^ have the behavior as n — ?■ oo 

Cn^-cn'', Dnr^dnP (39) 

with 

c=--(2-A2), /i = -l, rf=l, p = -2. (40) 

9 

By analysis similar to the 2-photon case, we have 

lim sup f|Q,T|n!)~ =0 (41) 

for all non-trivial solutions of (37). The zero limit means that the n > 1 part of the 
recurrence (16) possesses a minimal solution Q™", and the corresponding infinite power 
series expansions 0™"(z) generated by substituting Q™" for the Q^'s in (35) and (34), 
converge in the whole complex plane. 

The ratio of successive elements of the minimal solution Q™" can be expressed as 
continued fractions, 

^ Q™" _ At+l -Pn+2 -Pn+3 / .^^x 

"^ ~ Qmin ~ r r r -'" ' ^ ' 

which for n = reduces to 

5=^ = —^^^-^... r43) 

Wio ^\~ ^2~ ^Z~ 

0min 

The ratio 5*0 = Jlni„ involves QJT*", although the above continued fraction expression is 

rymin 

obtained from the 2nd equation of (37). On the other hand, the ratio 5*0 = nlmn of the 



first two terms of a minimal solution is unambiguously fixed by the n = part of the 
recurrence (37), that is, 



5*0 — —Co 



4gK, 



2kuj(2 -A^)-\E + u- A 



(44) 



In general, (43) and (44) can not be both satisfied for arbitrary values of the recurrence 
coefficients C„ and Dn- Thus general solutions to the recurrence (37) are dominant and 
are usually generated by simple forward recursion from a given value of Qq . Physical 
meaningful solutions are those that are entire in the Bargmann-Hilbert space [6]. They 
can be obtained if E can be adjusted so that equations (43) and (44) are both satisfied. 
Then the resulting solution sequence Q~{E) will be purely minimal and the power series 
expansion (34) will converge in the whole complex plane. Equating the right hand sides 
of (43) and (44) yields an implicit continued fraction equation for the regular spectrum 
E. That is, the regular energies E of the two- mode Rabi model are determined by the 
zeros of the transcendental function G{E) = Sq + Cq with 5*0 and Cq given by (43) and 
(44), respectively. Only for the denumerable infinite values of E which are the roots 
of G{E) = 0, do we get entire solutions of the differential equations (32) and (33). The 
transcendental eigenvalue equation G{E) = may be solved for E by standard root-search 
algorithms (see e.g. [21, 22] and references therein). 
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A Driven Rabi model 

Like the Rabi model, both the 2-photon and two-mode Rabi models studied in this paper 
have Z2 symmetry. However, our method can be applied to more general class of dynam- 
ical systems. As an example, in this appendix we revisit the driven quantum Rabi model 
[7] and derive a transcendental equation for its energy spectrum. 

The driven Rabi model is a generalization of the Rabi model with broken Z2 symmetry. 
It has recently been used in [23] to examine quantum thermalization. The Hamiltonian 
of the model reads 

HdR = ujb^b + Aa, + ga^{b^ + b) + 5a^, (45) 

where S is the drive amplitude. After a canonical Bogoliubov transformation b = a + 
A, b^ = a'^ + A, where A is a real parameter, the Hamiltonian (45) becomes 

HdR = cub^b + Aa, + {ga^ + \){b^ + b) + 6 a^ + 2g\ a^ + A=^. (46) 

Using the Bargmann realization a^ = z, a = ^, working in a representation defined by 
(Tx diagonal and choosing A = —g, we can turn the transformed driven Rabi Hamiltonian 

9 



into the matrix differential operator 



H, 



^z£ + 5-g' 



dR 



OJ \z 



dz 



2gz-5 + 3g^ 



(47) 



In terms of two-component wavefunction ip{z) = ((/9+(z), ip^{z))'^, the time-independent 
Schrodinger equation, HdR(p{z) = Eip{z), yields a system of couple of differential equa- 
tions 



u z—- + (5 -/-£'(/?+ + A V?- = 0, 



dz 



u z 



2g\ d 



CO J dz 



2gz + 3g'^ -5-E 



(/9_ + A(/9^ 



0. 



(48) 
(49) 



Solutions to these equations must be analytic in the whole complex plane if E belongs to 
the spectrum of Hf>. So we are looking for solutions of the form 



n=0 n=0 

which converge in the entire complex plane, i.e. solutions which are entire. 
Substituting (50) into (48), we obtain 

Thus IZ'^ is not analytic in E but has simple poles at 

-,2 



(50) 



E = nuj + 6 



9 



u 



n 



0,1,---. 



(51) 



(52) 



The energies (52) appear for special value of model parameters and correspond to the 
exceptional solutions of the driven Rabi model. If (52) is satisfied, the infinite series 
expansions (50) truncate and reduce to polynomials in z but only if the system parameters 
satisfy certain constraints. This can be easily verified by following the procedure in [19] 
(and thus the driven Rabi model is quasi-exactly solvable). Majority part of the spectrum 
of the driven Rabi model is regular for which (52) is not satisfied. The regular spectrum 
of the model is given by the zeros of the transcendental function Q{E) obtained below. 
From (49), we obtain the 3-term recurrence relation for TZ~ 



7^^ + Xo7^^ = o, 



0, n>l. 



(53) 



where 



X„ 



K 



1 



2^(n+l) 

1 
n + 1' 



3g2 

E-nuj + 5 —^ 

u E 



nu -6 + ^ 



(54) 
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The characteristic equation for the n > 1 part of the recurrence relation (53) ist'^ — j-t = 0, 
which gives two distinct roots ti = and ^2 = j-- Applying the Poincare-Perron theorem 
(i.e. theorems 2.1 and 2.2 of [20]), the two linearly independent solutions 7ln,i ^^^ ^n,2 
of the 2nd equation of (53) satisfy 

lim ^^ = t„ r = l,2, (55) 

n-i-oo 7?- r, , , V / 

and T?.™" = TZni is a minimal solution. 

The ratio of successive elements of the minimal solution sequence 7?.™" is expressible 
in terms of continued fractions, 

ipmin V V V 

rp _ l^n+l _ -'n+1 -'ra+2 -'n+3 /p-^x 

'<-n ^n+1— -^n+2— ^n+3^ 

which for n = reduces to 

^^WT^ ^J^J^ ^57^ 

'Tymin Y Y V ' ^ ' 

/Cq Ai— A2— A3 — 

However, for the single-ended sequences appearing in the infinite series expansions (50), 
the first equation (i.e. the n = part) of the recurrence (53) requires that 

T, = -X, = -^iE + 5-^-^+ f y (58) 

2g y u E-5 + ^-j 

In general, (57) and (58) can not be both satisfied for arbitrary values of the recurrence 
coefficients X„ and Yn- Physical meaningful solutions are those that are entire in the 
Bargmann-Hilbert space [6]. They can be obtained if E can be adjusted so that equations 
(57) and (58) are both satisfied. Then the resulting solution sequence TZn{E) will be 
purely minimal and the power series expansions (50) will converge in the whole complex 
plane. Equating the right hand sides of (57) and (58) yields a transcendental equation for 
the spectrum E. That is, regular energies of the driven Rabi model are given by the zeros 
of the transcendental function Q{E) = Tq + Xq, where Tq is the continued fraction in (57) 
and Xq is the right hand side of (58). The transcendental eigenvalue equation Q{E) = 
may be solved for E by standard root-search techniques. 
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